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Motivation (1/2) 


- Future aircraft systems will rely more on electrical and electronic 
components 
* UAV’s with all electric powertrain are increasingly being used for long 
missions 
* Electrical and Electronic components have increasingly critical role in 
on-board, autonomous functions for 
— Vehicle controls, communications, navigation, radar systems 


— Power electronic devices such as power MOSFETs and IGBTs are 
frequently used in high-power switching circuits 


— Batteries are the sole energy storage 
— The integrated navigation (INAV) module combines output of the GPS 
model and inertial measurement unit. 
* Assumption of new functionality increases number of faults with 
perhaps unanticipated fault modes 
* We need understanding of behavior of deteriorated components to 
develop capability to anticipate failures/predict remaining RUL 
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- Goals 
— Understand battery behavior through dynamic models 


— Develop model-based algorithms for state estimation, end of 
discharge (EOD) prediction, and end of life (EOL) prediction 


— Validate algorithms in the lab and fielded applications 
¢ Algorithms 

— Prognostic Architecture 

— Dynamic state and state-of-charge estimation 
¢ Modeling 

— Electric circuit equivalent (for EOD prediction) 

— Electrochemistry-based model (for EOD and EOL prediction) 
* Applications 

— Rover 

— Edge 540-T electric aircraft 
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Outline 


¢ Introduction to Prognostics 

¢ Introduction to Model-based Prognostics 
¢- Research Approach 

¢ Electric Vehicle Powertrain 


¢ Algorithms 

— Prognostic Architecture 

— Dynamic state and state-of-charge estimation 
* Applications 

— Rover 

— Edge 540-T electric aircraft 
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Integrated Prognostics Architecture 


- System (battery) gets inputs (current) and produces outputs (voltage) 


- State estimation computes estimate of state given estimates of age 
parameters 


¢- EOD prediction computes prediction of time of EOD, given state and 
age parameter estimates 


- Age parameter estimation computes estimates of age parameters 

- Age rate parameter estimation computes parameters defining aging 
rate progression 

¢- EOL prediction computes prediction of time of EOL, given age 
parameter and age rate parameter estimates 
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The Basic 


Threshold as a Function of System State 
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Why Prognostics? 


Example: UAV Mission 


Visit waypoints to accomplish science objectives. Predict aircraft battery end of discharge to 


determine which objectives can be met. Based on prediction, plan optimal route. Replan if 
prediction changes. 
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Why Prognostics? 


¢ Prognostics can enable: 


Adopting condition-based maintenance strategies, instead of time- 
based maintenance 


Optimally scheduling maintenance 
Optimally planning for spare components 


Reconfiguring the system to avoid using the component before it 
fails 


Prolonging component life by modifying how the component is used 
(e.g., load shedding) 


Optimally plan or replan a mission 


- System operations can be optimized in a variety of ways 
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The Basic Idea : Batteries Example 
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The Basic Idea : Batteries Example 
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———Prognostic Algorithm Categories 


. TYPE I: Reliability Data-based 
Use population based statistical model 
— These methods consider historical time to failure data which are used to model 
the failure distribution. They estimate the life of a typical component under 
nominal usage conditions. 


— Ex: Weibull Analysis 


. wpe ll: Stress-based 
Use population based fault growth model — learned from accumulated knowledge 
— These methods also consider the environmental stresses (temperature, load, 
vibration, etc.) on the component. They estimate the life of an average 
component under specific usage conditions. 


— Ex: Proportional Hazards Model 


¢ Type Ill: Condition-based 
— Individual component based data-driven model 
— These methods also consider the measured or inferred component degradation. 
They estimate the life of a specific component under specific usage and 
degradation conditions. 
— Ex: Cumulative Damage Model, Filtering and State Estimation 


Physics-Based Methods 


¢ Description of a system’s underlying physics using suitable 
representation 
¢ Some examples: 
— Model derived from “First Principles” 
* Encapsulate fundamental laws of physics 
= PDEs 
= Euler-Lagrange Equations 


— Empirical model chosen based on an understanding of the dynamics of a 
system 


¢ Lumped Parameter Model 

¢ Classical 15 (or higher) order response curves 
— Mappings of stressors onto damage accumulation 

¢ Finite Element Model 

¢ High-fidelity Simulation Model 


¢ Something in the model correlates to the failure mode(s) of interest 


Physics-Based Models 


¢ Pros 
— Results tend to be intuitive 
¢ Based on modeled phenomenon 
¢ And when they’re not, they’re still instructive (e.g., identifying needs for more 
fidelity or unmodeled effects) 


— Models can be reused 
¢ Tuning of parameters can be used to account for differences in design 


— If incorporated early enough in the design process, can drive sensor 
requirements (adding or removing) 


— Computationally efficient to implement 


¢ Cons 
— Model development requires a thorough understanding of the system 
— High-fidelity models can be computationally intensive 


¢ Examples 
— Paris-Erdogan Crack Growth Model 
— Taylor tool wear model 
— Corrosion model 
— Abrasion model 


INTRODUCTION TO MODEL- 
BASED PROGNOSTICS 


Model-based prognostics (1/2) 


_ ¢ State vector includes 
a) 7 Le), u(t) ‘ my) dynamics of the degradation 
y(t) = h(x(t)), u(t)) + v(k) process 


It might include nominal 
operation dynamics 
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Model-based prognostics (2/2) 


¢ Tracking of health 
state based on 
measurements 


- Forecasting of health 
state until failure 
threshold is crossed 

- Compute RUL as 
function of EOL 
defined at time failure 
threshold is crossed 
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RESEARCH APPROACH 


High level research efforts 


- Prognostics models and algorithms 
— Identification of precursors of failure for MOSFETs under different failure mechanism 
conditions 
— ldentification of precursors of failure for different IGBT technologies (CALCE) 
— Modeling of degradation process MOSFETs 
— Development of prognostics algorithms 
- Prognostics for output capacitor in power supplies (Vanderbilt University) 
— Electrical overstress and thermal overstress 
— Development of prognostics algorithms 
- Accelerated Life Testing 
— Thermal overstress aging of MOSFETs and IGBTs 
— Electrical overstress aging testbed MOSFETs 
— Electrical overstress aging testbed for Capacitors 
- Effects of lightning events of MOSFETS (LaRC) 
- Battery Degradation and ageing ( ARC — LaRC) 
- Ageing Effecting on ESC’s ( ARC — LaRC) 
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Research Approach 


Kol=Yahihixersid(o)ame)mr-lilUlcomanveye(stswr-lalemiarsyima)t-lice)atsjal/emcomsals)/mey-latceult-lanrclilelas) 
laatsxeiarelalisiaats 


Ko f=Yalihiter-\ile)ale) ime) =\e1 ese) ese) mre lilU leo m'2al (evalu o)t-WVar-lam>ssss\>)a0it- lM ce)(om lam tate e)g=vol(eil(e)ame)inc=yant-lialialemilic 


IDY=\V=)(o) 0) aat=1al me) m= (eex>)(=1e-1¢2Ye m= (e]| ale mm (ossjhel-re lcm ial-lmr-lelllitclomlal=m=>.<0)(e)e-ile)ame)melliicic-lalmr-Ullll a= 
faat=vevat-lalisjantsw-larem-l(omiat>m0lale(s1esit-larel|aleme)imelt=lpar-\el>me)cele]g=sss-1(0)a) 


IDY=W{=1(o) ola al= al meme (=e |e-(er-tile)amanrele(=)\om er-ks\-Le me) alm alm e)ahvss) (eco) miat=me(-\Uleroi- [ale mialo mri Ula 
laatsxeiarslalisiaats 


IDY=WV/=)(o) o)pat=)a) Movin c= Janie Uialialem li=m e)a=xe|(eqt(e)am-l(e(e)aiialaatcmuarumr-l.¢om al (om-(ecerelUlalmial-Melli(c)t-almcyelU|gerc\oMe) | 
Ui ater=}at-lialavadaliiomcs\ds1e-(e]|ale mm elaniAsiiercss ef=ts\o1e me (cle) e-(er-1i(e)amanlele(=)\omlar-imexe)alci(e(=]ecmie lie | come) el=l¢-1i(e) atl 
reValom=)anviice)alaaisyaltclimexe)ale lise) ats 


Prognostics Center of Excellence 


Prognostics Algorithm Maturation through Validation 
Experiments 
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Prognostics Algorithm Maturation through Validation 
Experiments 
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ELECTRIC VEHICLE 
POWERTRAIN 


System Level Prognostics 
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Hardware in Test Loop 
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Accelerated Aging 


- Traditionally used to assess the reliability of products with expected 
lifetimes in the order of thousands of hours 
— in aconsiderably shorter amount of time 
- Provides opportunities for the development and validation of prognostic 
algorithms 
* Such experiments are invaluable since run-to-failure data for 
prognostics is rarely or never available 
* Unlike reliability studies, prognostics is concerned not only with time to 
failure of devices but with the degradation process leading to an 
irreversible failure 
— This requires in-situ measurements of key output variables and observable 
parameters in the accelerated aging process with the associated time 
information 
¢ Thermal, electrical and mechanical overstresses are commonly used 
for accelerated aging tests of electronics 
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State Estimation 


¢ What is the current system state and its associated 
uncertainty? 
— Input: system outputs y from k, to k, v(k,:k) 
— Output: p@x(k), AK)|Wko:k)) 

¢ Battery models are nonlinear, so require nonlinear state 
estimator (e.g., extended Kalman filter, particle filter, 
unscented Kalman filter) 


- Use unscented Kalman filter (UKF) 
— Straight forward to implement and tune performance 


— Computationally efficient (number of samples linear in size of state 
space) 
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Prediction 


- Most algorithms operate by simulating samples forward in 
time until £ 
- Algorithms must account for several sources of uncertainty 
besides that in the initial state 
— A representation of that uncertainty is required for the selected 
prediction algorithm 
— A specific description of that uncertainty is required (e.g., mean, 
variance) 
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Prediction Algorithm 


* The P function takes an initial state, — Ajgorithm 1 kpikp) — Pixikp).O,,.Uz,.Vn.) @x,.Uxp. Vip) 
and a parameter, an input, and a 
process noise trajectory 

— Simulates state forward using f until £ is 
reached to computes 4, for a single 
sample 

* Top-level prediction algorithm calls P 

— These algorithms differ by how they 
compute samples upon which to call P 


- Monte Carlo algorithm (Mc) takes as 


ke kp 

x(k) — x(kp) 

while Tr (x(k), Ox%,(k), Ur, (k)) = Odo 
x(k +1) + £(k, x(k), Ox, (k), Urp(k), Vep(F)) 
kek+1 
x(k) — x(k +1) 

end while 

Kkpe(kp) ek 
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Input Algorithm 2 {kj}, = MC(p(x(kp), @(kp)|y(ko:kp)), 
— Initial state-parameter estimate P(A@). pP(Au). P(Av) NV) 
— Probability distributions for the i: fori =1to N do 
surrogate variables for the parameter, 2 (x(kp),A(kp)) ~ p(x(kp), 0(kp)|y(ko:kp)) 
input, and process noise trajectories 3 AO w p(Ag) 


— Number of samples, V 


¢ McC samples from its input 
distributions, and computes k,, 


3 
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5: Au ~ p(Au) 
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POWER TRANSISTORS 


Modeling for Power MOSFET under electrical overstress 


¢ Two-transistor model is shown to be a 
good candidate for a degradation model 
for model-based prognostics. 

¢ The model parameters K, and W1 could 
be varied as the device degrades as a 
function of usage time, loading and 
environmental conditions. 


- Parameter W1 defines the area of the a 
healthy transistors, the lower this area, the - 
larger the degradation in the two-transistor 
model. In addition, parameter K serves as Nig, ree ek 


a scaling factor for the thermal resistance 
of the degraded transistors, the larger this 
factor, the larger the degradation in the 
model. 
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Precursor of Failure 


- As case temperature increases, ON- 
resistance increases 

- This relationship shifts as the 
degradation of the device increases 

- For adegraded state, ON-resistance 
will be higher at any given case 
temperature 

- This is consistent with the die-attach 
damage since it results on increased 
junction temperature operation 


¢ This plot can be used directly for 


R vs. Temperature: Device 9 
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Degradation process data 


Normalized ON-state resistance (ARpgon)) 
and filtered trajectory for device #36 


— RDS,,, 
Filtered 


ARDS, 


Cases #08, #09, #11, #12 and #14 
are used for algorithm development 
purposes. 

Case #36 is used to test the 
algorithms. 


Normalized ON-state resistance (ARpg(on)) 
and filtered trajectory for device #36 
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Empirical Degradation Model 


¢ An empirical degradation model was selected for the model- 
based algorithms 


¢ Exponential based function to capture degradation process 
* Two parameters in the model which will be estimated on-line 
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Prediction of Remaining Life 


RUL Prediction Methodology Considerations 


A single feature is used to assess the health state of the device 
(ARpscon)) 

It is assumed that the die-attached failure mechanism is the only active 
degradation during the accelerated aging experiment 

Furthermore, ARpgon accounts for the degradation progression from 
nominal condition through failure 

Periodic measurements with fixed sampling rate are available for 
ARpscon) 

A crisp failure threshold of 0.05 increase in ARpgion is used 

The prognostics algorithm will make a prediction of the remaining 
useful life at time ¢,, using all the measurements up to this point either 
to estimate the health state at time tf, in a regression framework or in a 
Bayesian state tracking framework 

It is also assumed that the future load conditions do not vary 
significantly from past load conditions 
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RUL Prediction Algorithms 


* Gaussian Process Regression 
— Algorithm development cases used to select covariance matrix structure 
and values 
¢ Extended Kalman filter 
— Empirical degradation model 


— State variable: Normalized ON-resistance and degradation model 
parameters 


— Arbitrary values for measurement and process noise variance 
* Particle filter 
Empirical degradation model 
State variable: Normalized ON-resistance, degradation model parameters 
Exponential growth model used for degradation model parameters 
Arbitrary values for measurement and process noise variance 
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RUL estimation results 
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Degradation Mechanisms 


Electrolyte Evaporation 
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Decrease in capacitance 
Increase in ESR 


Failure Modes 
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Degradation Model: Electrical Circuit Equivalent 
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Capacitance Degradation Model 


Decrease in electrolyte volume : 


(1) 
Ve(t) = Veo — (WeAsJeot) 
where: 
V: dispersed volume at time t, V.: initial electrolyte volume 
A,: surface area of evaporation, je, : evaporation rate 
t : time in minutes, w. = volume of ethyl glycol molecule 
Capacitance (C) ): Physics-Based Model: 
(2) 


C = (2ereEoAs)/dco 


Electrolyte evaporation dominant degradation phenomenon 
— First principles: Capacitance degradation as a function of electrolyte loss 


2€ R€0 V0 = V.(t) 3 
where: D, : C(t) = ( ———e ; 2 
C : capacitance of the capacitor, dc Jeo t We 


€p : relative dielectric constant, 


€Q : permittivity of free space, 
de : oxide thickness. 
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Capacitance Degradation Model 


* Oxide breakdown observed - experimental data 
- The breakdown factor is exp. function of electrolyte evaporation 


Cox) = EXP f(Veo — Vert) 


* Updated in capacitance degradation model : 


C= (2€ Rep AS CHK )/de, 


: 2€ REC Veo — V.(t) 
Di HC (t) = Cbk(t) ( x ) (a) 


Prognostics Center of Excellence 


44 


Dynamic Model of ESR 


¢ Decrease in electrolyte volume : 
Ve(t) = Veo — (weAsJeot) 
- ESR 
— Based on mechanical structure and electrochemistry. 
— With changes in R, (electrolyte resistance ) 
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RUL and Validation - EOS -Experiment — ESR 


Degradation Model D; 


Tracking 


Alpha Lambda 


RUL (hours) 


Cap #2 


—s— measured data 
+ filter data 


a 
re: 


L L L L L L 
80 100 120 140 160 180 
Aging Time ( Hours) 
Output Error - Cap #2 


t t : i t 7 t 
80 100 120 140 160 180 200 
Aging Time ( Hours) 


a=0.3, B=0.5 


[ c 
140 160 


t t 
180 200 
Le 149 


Predictions at different aging time 


: 
160 


[ [ 
180 200 
L.= 171 


[ c c [ c 
80 100 120 140 160 
Aging Time ( Hours) 


go 100 120 
Aging Time (Hours) 


Prognostics Center of Excellence 


[ [ 
180 200 
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LI- ION BATTERIES 


Electrochemical Li-ion Model 


= Electrochemical Models vs. Empirical Models 
Battery physics models enable more direct representation of age-related changes in 
battery dynamics than empirical models 
» Typically have a higher computational cost and more unknown parameters 


- Lumped-parameter, ordinary differential equations 


¢ Capture voltage contributions from different sources 
— Equilibrium potential >Nernst equation with Redlich-Kister 


expansion 

— Concentration overpotential > split electrodes into surface and bulk 
control volumes ; ss 

— Surface overpotential > : — Measured 
Butler-Volmer equation ‘ ~ = = Predicted 


applied at surface layers 


— Ohmic overpotential > 
Constant lumped resistance 
accounting for current - 
collector resistances, 324 | | | 
electrolyte resistance, 0 L000 2000 3000 on 5000 6000 7000 
solid-phase ohmic resistances |. -ocience 


ognostics 


Voltage (V) 


Battery Aging 


- Contributions from both decrease in mobile 4 — Simulated 
Li ions (lost due to side reactions related to = > 
aging) and increase in internal resistance i | " 


— Modeled with decrease in “g”*” parameter, 
used to compute mole fraction 

— Modeled with increase in “Rk,” parameter rr rr re 
capturing lumped resistances biaccee oe 


Cycle 16 
Cycle 26 
Cycle 36 


Cycle 46 : | \ —— 


i Measured 
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Cycle 76 
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~~ Cycle 146 
— Cycle 156 
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Rover 


- Planetary rover testbed at NASA Ames 
Research Center 


— 24 lithium ion batteries, two parallel sets of 
12 in series 


— Batteries power 4 motors, one for each 
wheel (skid steering) 

¢ Rover operated in two driving modes 
— Unstructured driving 


- Rover is driven freely by an operator, 

without prior knowledge of actions 
— Structured driving 

- Rover has a given mission, to visit a 
set of waypoints 

* Rover moves along, visiting waypoints 

- End-of-discharge prediction is 
required in order to ensure the given 
set of waypoints can be visited, and if 
not, to replan the route to optimize 
mission value 


Ref : A. Sweet et al “Demonstration of Prognostics-Enabled Decision Making Algorithms on a 
Hardware Mobile Robot Test Platform’, PHM 2013 
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Results: Structured Driving 


x Rover Path 
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Predictions are very accurate since rover travels at a known fixed average speed, 


and waypoints are known. 


Uncertainty in predictions is significantly less than for unstructured driving, since 
more information about future inputs are known. 
Predictions are under at the start because power drawn for first 500 s is half the 


average. 
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Edge 540-T 


Subscale electric 
aircraft operated 

at NASA Langley 
Research Center 


Powered by four 
sets of Li- 
polymer batteries 


Estimate SOC 
online and 
provide EOD and 
remaining flight 
time predictions 
for ground-based 
pilots 


é 
Throttle 


13 4 
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Predication over Flight Plan 


Predicted Time Remaining (min) 


¢ Measured and predicted 


¢ The min, max and median 
predictions are plotted 
from each sample time 
until the predicated SOC 
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¢ Predictions for remaining flight time for 
entire flight plan 

¢ Overestimate till parasitic load is injected 

¢ Once the parasitic load is detected the 
remaining flying time time prediction shifts 


5 10 


Time (min) 


= down. 
Ref : E. Hogge et al, “Verification of a Remaining Flying Time Prediction System for Small Electric Aircraft”, PHM 2015 


Prognostics Center of Excellence 


Performance Requirements 


¢ Accuracy requirements for the two minute warning were specified as: 

— The prognostic algorithm shall raise an alarm no later than two minutes 
before the lowest battery SOC estimate falls below 30% for at least 90% of 
verification trial runs. 

— The prognostic algorithm shall raise an alarm no earlier than three minutes 
before the lowest battery SOC estimate falls below 30% for at least 90% of 
verification trial runs. 

— Verification trial statistics must be computed using at least 20 experimental 
runs 
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Data Sets Available for Download 


- https://ti.arc.nasa.gov/tech/dash/pcoe/prognostic-data-repository/ 


Randomized Battery Usage Data Set 
Publications using this data set HIRF Battery Data Set 


Batieries are continuously cycled with randomly generated current profiles. | Publications using this data set 
Reference charging and discharging cycles are also performed after a fixed 
interval of randomized usage in order to provide reference benchmarks for 
battery state of health. 


Description Battery Data collected from the Experiments on the Edge 540 Aircraft in HIRF 
(Chamber. Refernce document can be downloded here 
Format The set is in .mat format and has been zipped. 
Datasets + Download HIRF Battery Data Set 1 (184 downloads) 
+ Download HIRF Battery Data Set 2 (127 downloads) 
+ Download HIRF Battery Data Set 3 (131 downloads) 
+ Download HIRF Battery Data Set 4 (125 downloads) 
+ Download HIRF Battery Data Set 5 (149 downloads) 
+ Download HIRF Battery Data Set 6 (135 downloads) 
Dataset _C. Kulkarni, E. Hogge, C. Quach and K. Goebel "HIRF Battery Data Set", 
Citation NASA Ames Prognostics Data Repository 
(http://ti.arc.nasa.gov/project/prognostic-data-repository), NASA Ames 
|Research Center, Moffett Field, CA 
Publication Verification of a Remaining Flying Time Prediction System for Small Electric 
Citation Aircraft. Edward F. Hogge, Brian M. Bole, Sixto L. Vazquez, Jose R., Annual 
Conference of the Prognostics and Health Management, PHM 2015 


le Dowslons Reacontaed Satery Usage Oates Set 1 (1285 downloads) 


at 2 (936 downloads) 
+ Download Randomized Battery Usage Data Set 3 (906 downloads) 
+ Download Randomized Battery Usage Data Set 4 (4217 downloads) 
+ Download Randomized Battery Usage Data Set 5 (825 downloads) 
+ Download Randomized Battery Usage Data Set 6 (890 downloads) 
+ Download Randomized Battery Usage Data Set 7 (857 downloads) 
B. Bole, C. Kulkarni, and M. Daigle "Randomized Battery Usage Data Set", 
NASA Ames Prognostics Data Repository 
(httpy/ti.arc.nasa.gov/project/prognostic-data-repository), NASA Ames 
Research Center, Moffett Field, CA 
B. Bole, C. Kulkarni, and M. Daigle, ‘Adaptation of an Electrochemistry-based 
Li-lon Battery Model to Account for Deterioration Observed Under 
Randomized Use’, Annual Conference of the Prognostics and Health 
Management Society, 2014 
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Remarks (1/2) 


¢ Electrical and Electronics PHM Maturity - scientific 
and engineering challenges 


¢« Research approach challenges 


— How to balance lack of knowledge of the system vs own 
expertise on particular PHM tools 


— Data-driven or model-based? 

- Data is always needed but more important, 
information about degradation/aging processes 
is key 

¢ Experiments and field data 
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Remarks (2/2) 


- Aging systems as a research tool 


— Value in terms of exploration of precursors of failure and 
their measurements is evident 


— Still an open question on how degradation models and 
algorithms are translated to the real usage timescale 


- In the use of physics 
— It should be embraced 

- Validate models and algorithms with data from lab 
experiments and fielded systems 


- A success in developing PHM methodologies in an 
real usage application will require the right team 
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